Rotational energy transfer in H2+H2 
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Quantum-mechanical close-coupling calculations for state-to-state cross sections and thermal rates are reported 
for H2+H2 collisions. Two recently developed potential energy surfaces (PES) for the H2— H2 system are applied, 
namely, the global potential surface from the work of A.I. Boothroyd, P.G. Martin, W.J. Keogh, M.J. Peterson, 
J. Chem. Phys., 116 (2002) 666, and a restricted, model surface from the works of P. Diep, J.K. Johnson, J. 
Chem. Phys., 113 (2000) 3480; ibid. 112, 4465. The low temperature limit is investigated. We found significant 
differences in cross sections and corresponding thermal rates calculated with these two PESs. 



1. Introduction 

The investigation of elastic and inelastic colli- 
sions between molecules and between molecules 
and atoms can provide valuable information 
about interactions, chemical properties and en- 
ergy transfer dynamics [1-15]. The hydrogen 
molecule is the simplest and most abundant 
molecule in the universe's molecular clouds and 
plays an important role in many areas of as- 
trophysics. For example, knowledge of the ro- 
vibrational excitation and de-excitation rate con- 
stants in molecular hydrogen collisions is of 
fundamental importance for understanding and 
modeling the energy balance in the interstellar 
medium. The energy transfer processes involv- 
ing H2 molecules control the evolution of shock 
fronts and photodissociation regions (PDRs) in 
the interstellar medium. Additionally, the en- 
ergy transfer between H2 molecules and between 
H2 and other atoms/molecules is important for 
cooling of primordial gas and shock wave-induced 
heating in the interstellar media. However, to 
accurately model the thermal balance and kinet- 
ics of such important systems one needs accurate 
statc-to-state rate constants k 



VJV'J 



(T). 



Experimental measurement of quantum state 
resolved cross sections and rates is a very difficult 
technical problem. On the other hand, accurate 
theoretical data requires precise potential energy 



surfaces and reliable dynamical treatment of the 
collision processes. The first attempt to construct 
a realistic full-dimensional ab initio PES for the 
H2— H2 system was done in works [2], and the 
potential was widely used in a variaty of methods 
and computation techniques. 

Recently the H2— H2 system has been reinvesti- 
gated and an accurate interaction potential from 
the first principles has been developed in work 
[12| . However, in this work the Diep and Johnson 
potential energy surface (DJ PES) was extrapo- 
lated only for the rigid rotor monomer model of 
H2— H2. 

An extensive study of the H2— H2 PES has been 
recently reported by Boothroyd et al. |13| . where 
the potential energies have been represented at 
48180 geometries respectively with a large basis 
set at the multireference configuration interaction 
level. 

In this work we provide a test of sensitivity for 
the new PESs for collisions of rotationally excited 
H 2 molecules 



H 2 (ii)+H 2 (i 2 )^H 2 (j 1 )+H 2 (i 2 ; 



(1) 
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We apply only the new surfaces mentioned above. 

The scattering cross sections and their corre- 
sponding rate coefficients are calculated using a 
non reactive quantum-mechanical close-coupling 
approach. In the next section we will briefly out- 
line the method. Our results and discussion are 
presented in Section 3. 
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2. Method 

In this section we provide a brief outline of 
the quantum-mechanical close-coupling approach 
used in our calculations. All relevant details 
have been developed in works £Q. The hydrogen 
molecules are treated as linear rigid rotors. The 
model has been applied in few previous works, see 
for example |1I4| . For the considered range of ki- 
netic energies of astrophysical interest the rotor 
model is considered to be adequate |T3]. 

The cross sections for rotational excitation and 
relaxation phenomena can be obtained directly 
from the S'-matrix. In particular the cross sec- 
tions for excitation from j\j 2 3x3% summed 
over the final m!-yVn! 2 and averaged over the ini- 
tial mim2 corresponding projections of the H2 
molecules angular momenta j\ and j 2 are given 
by 

<r{3x,3%,3x3%,e) - (2il + 1)(2i2 + 1)Iw 

(2J+l)\S aa ,-SL,(E)\ 2 . (2) 

Jji2j' 12 LL' 

The kinetic energy is e = E — Bxjx(jx + 1) — 
B 2 j 2 (j 2 + !)• Here E is the total energy in 
the system, Bu 2 \ = 60.8 cm -1 are the rota- 
tion constants of the colliding H2 molecules, J 
is total angular momenta of the 4-atomic sys- 
tem, a = (jx323x%L), where j x + j 2 = jx% and 
J12 + L = J, k aa , = 2M 12 (E + E a - E a ,)V 2 
is the channel wavenumber and E a r a i\ are rota- 
tional channel energies. 

The log-derivative matrix is propagated to 
large i?-intermolecular distances, since all exper- 
imentally observable quantum information about 
the collision is contained in the asymptotic be- 
haviour of functions U^ M (R — > 00). A set of cou- 
pled second order differential equations for the 
unknown radial functions U^ M (R) is solved 

( d 2 L(L + l) l2 \ TT j M(T} s 0M 
E / <<l>t M (h,h,R)\V(r l ,r 2 ,R)\ 

a' 

^(h,f 2 , R) > UiF{R)df 1 df 2 dR, (3) 
where a = {j\j 2 j\ 2 L). We apply the hybrid mod- 



ified log-derivative- Airy propagator in the general 
purpose scattering program MOLSCAT ^H] to 
solve the coupled radial equations J2J- Different 
propagator schemes included in MOLSCAT have 
been tested. Our calculations showed that other 
propagators can also provide quite stable results. 

The numerical results are matched to the 
known asymptotic solution to derive the physi- 
cal scattering S'-matrix 

T/J A , p -Kk aa R-(ln/2)) _ 

( h^2L \ 1 QJ p -i{k aa ,R-{!' ■*/%)) (a\ 
\k aa ,J * aa ' e ' 14 > 

The method was used for each partial wave until a 
converged cross section was obtained. It was veri- 
fied that the results are converged with respect to 
the number of partial waves as well as the match- 
ing radius, R m ax, for all channels included in our 
calculations. 

The relationship between the rate coefficient 
kj 1 j 2 ^j' i j' 2 (T) and the corresponding cross section 
a jih^i' 1 j' 2 (Ekin) can be obtained through the fol- 
lowing weighted average 

, frp . _ 8k B T 1 [°° 

Khh^m^ ) - — {kBT)2 

°h^j' 1 3 2 (e)e- £/kBT , (5) 

where e = E — Ej 1 — Ej 2 is precollisional transla- 
tional energy at the translational temperature T, 
ks is Boltzman constant, /x is reduced mass of the 
molecule-molecule system and e s is the minimum 
kinetic energy for the levels ji and j 2 to become 
accessible. 

3. Results 

As we mentioned in the Introduction, in this 
work we apply the new PESs from works |12j and 
[T3| . The DJ PES Q2j is constructed for the vi- 
brationally averaged rigid monomer model of the 
H 2 — H 2 system to the complete basis set limit us- 
ing coupled-cluster theory with single, double and 
triple excitations. A four term spherical harmon- 
ics expansion model was chosen to fit the surface. 
It was demonstrated, that the calculated PES can 
reproduce the quadrupole moment to within 0.58 
% and the experimental well depth to within 1 %. 
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The bond length was fixed at 1.449 a.u. or 
0.7668 A. DJ PES is defined by the center-of-mass 
intermolecular distance, R, and three angles: Q\ 
and 82 are the plane angles and 4>i2 is the relative 
torsional angle. The angular increment for each of 
the three angles defining the relative orientation 
of the dimers was chosen to be 30°. 

The BMKP PES [HI is a global six-dimensional 
potential energy surface for two hydrogen 
molecules. It was especially constructed to repre- 
sent the whole interaction region of the chemical 
reaction dynamics of the four-atomic system and 
to provide an accurate as possible van der Waals 
well. In the six-dimensional conformation space 
of the four atomic system the conical intersection 
forms a complicated three-dimensional hypersur- 
face. Because the BMKP PES uses cartesian 
coordinates to compute distances between four 
atoms, we have devised some fortran code, which 
converts spherical coordinates used in Sec. 2 to 
the corresponding cartesian coordinates and com- 
putes the distances between the four atoms. In 
all our calculations with this potential the bond 
length was fixed at 1.449 a.u. or 0.7668 A as in 
the DJ PES. 

A large number of test calculations have also 
been done to secure the convergence of the re- 
sults with respect to all parameters that enter 
into the propagation of the Schrodinger equation 
[TT] . This includes the intermolecular distance 
R, the total angular momentum J of the four 
atomic system, Ni v i the number of rotational lev- 
els to be included in the close coupling expan- 
sion and others (see the MOLSCAT manual |15|h 
We reached convergence for the integral cross sec- 
tions, <r(Ekin), in all considered collisions. In the 
case of the D J PES the propagation has been done 
from 2 A to 10 A, since this potential is defined 
only for those specific distances. For the BMKP 
PES we used r min = 1 A to r max = 30 A. We 
also applied a few different propagators included 
in the MOLSCAT program. 

Now we present our results for differ- 
ent rotational transitions in collisions between 
para/ para- and ortho-/ortho- hydrogen molecules 
JI). Our main goal in this work is first: to 
carry out complete quantum-mechanical calcula- 
tions with new potentials for different transitions 
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Figure 1. Rotational state resolved integral cross 
sections for elastic scattering in the case of para- 
/para- and ortho-/ortho-hydrogen and transi- 
tions, when ji = ]2 = — ► j[ = 2, j 2 = 2 and 
ji = j2 = 1 — > j[ = 1, j 2 = 3- Calculations are 
done with the DJ (bold lines) and BMKP (open 
triangles up) PESs. The compensating factor of 
2 is included only in the elastic cross sections. 



in P-H2-I-P-H2 and 0-H2+0-H2 collisions and, sec- 
ond, to provide a comparative study and check 
the sensitivity of the two newest surfaces for the 
very important and fundamental H2-H2 system. 

The energy dependence of the elastic integral 
cross sections on the total energy (j e /(_Z?^ ^) in 
the H2+H2 system is represented in Fig. 1 (up- 
per plots) together with the state-resolved inte- 
gral cross sections Cj 1 j 2 ->j[j' 2 (Etot) for the ji = 
h = -> j[ = 2,j 2 = 2 and ji = j 2 = 1 -> 
j[ = l,j' 2 = 3 rotational transitions (lower plots) 
for both the BMKP and the DJ PESs respec- 
tively. As can be seen both PESs provide the 
same type of the behaviour in the cross section. 
These cross sections are in basic agreement with 
recent time-dependent quantum-mechanical cal- 
culations in work 7 . However, our results show, 
that the DJ PES generates about 30 % higher val- 
ues for the cross sections relatively to the BMKP 
PES. 

Also, it is important to point out here, that for 
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comparison purposes we don't include the com- 
pensating factor of 2 mentioned in 2 . However, 
in Fig. 2 (left plot) and in our subsequent calcula- 
tions of the thermal rate coefficients, kjj>(T), the 
factor is included. 

Significant differences in the cross sections of 
the two potentials are reflected in the state- 
resolved transition states j\ = 0, ji = — > ,][ = 
2, j' 2 = 0, as shown in Fig. 2 (right panel). That is 
why it seems that the DJ PES can provide much 
better results, as seen in the same figure in the left 
panel. Specifically, when we present the results 
for the corresponding thermal rates koo-2o{T) 
calculated with the DJ potential together with 
results of other theoretical calculations 0] and 
[5] the agreement is almost perfect. Next, Fig. 
3 provides the same results for the cross sections, 
but includes the details at low energies, together 
with the corresponding experimental data from 
U3- Now one can better see the considerable dif- 
ferences in these two cross sections. However, the 
DJ PES is able to provide very good agreement 
with experimental data Thus we conclude, 
that DJ PES is much better suited for the H 2 — H 2 
system. Moreover, in Fig. 4 we provide thermal 
rates for different transition states calculated with 
only the DJ PES and in comparison with other 
theoretical data obtained within different dynam- 
ical methods and PESs. Again the agreement is 
very good. 

Finaly, Table 1 depicts the thermal rates 
koo^2o{T) and kg2^22(T) at lower kinetic tem- 
peratures, up to 10 K. Numbers in each column 
should be multiplied by corresponding prefactor 
P. As can be seen our results are in very good 
agreement with recent experimental and theoret- 
ical data at higher temperatures. However, at 
T < 30 K we could not reproduce effectively the 
results of [HI ■ It may be because of the interatomic 
(H2-H2) distance restrictions in the DJ surface in 
which the potential ranges from 2 A to only 10 
A. It seems plausible, that for such low energies 
one needs to carry out calculations even beyond 
the 10 A limit. 

We provide close-coupling quantum- 
mechanical calculations of the state-resolved 
rotational excitation and deexcitation cross sec- 
tions and rates for molecular hydrogen collisions. 




Figure 2. Temperature dependence of the state- 
resolved thermal rate constant (left panel) and 
corresponding cross sections (right panel) for the 
transition j\ = j% = — * j[ = 2,j' 2 = 0. Squars 
and circles are the data from |3] and |S] respec- 
tively. Our results with the DJ PES are depicted 
with solid lines, the open triangles up represent 
our results with the BMKP surface. 
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• Experiment: B. Mate, et al., JCP 122, 0643 13 (2005) 
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Figure 3. Cross sections for the 00— >20 rotational 
transition calculated with the DJ and BMKP 
PESs for the H2+H2 collision. The circles are 
some experimental data from the work , trian- 
gles up and triangles down are the results of this 
work using the DJ and BMKP PESs respectively. 



Rotational energy transfer 



5 



A test of convergence and the results for cross 
sections and rate coefficients using two different 
potential energy surfaces for the H 2 — H 2 system 
have been obtained for a wide range of kinetic 
energies [TT] . 

Our calculations revealed, that both PESs can 
provide the same type of behaviour in regard 
to cross sections and rates for different transi- 
tion states. However, significant differences in 
the cross section of the 00— > 20 transition have 
been found. Also, it was already indicated, that 
at even larger kinetic energies the DJ potential 
overestimates relative to the BMKP surface the 
cross sections and rates by about 20-40 % 

In conclusion, the results of these calculations 
show, that additional work is needed to further 
improve the BMKP PES, particularly the part 
of the surface, which is responsible for the 00— ► 
20 transition. We also found, as can be verified 
from Table 1, it is difficult to use the DJ surface 
and carry out reliable calculations at tempera- 
tures lower than 30 K. 



Table 1 

Thermal rates fcoo^2o(T) & n d fco2^22(? 1 ) (m 3 s _1 ) 
at low and very low temperatures calculated 
with the DJ surface in comparison with available 
experimental*^ and theoretical data from |5] 







^00^20 




k 2 


-^22 




T(K) 


DJ 


9 * 


M 


DJ 




P 


100 


1.64 


2.2(4) 


1.92 


2.63 


3.23 


10 -au 


60 


3.81 


6.0(7) 


4.75 


5.62 


7.97 


10 -22 


30 


0.54 


1.1(1) 


0.80 


0.65 


1.34 


10 -2S 


20 


0.87 


2.3(3) 


1.63 


0.93 


2.71 


1Q -29 


10 


0.38 


2.7(7) 


1.65 


0.33 


2.62 


10 -40 



REFERENCES 

1. S. Green, J. Chem. Phys., 62 (1975) 2271. 

2. G. Danby, D.R. Flower, T.S. Monteiro, Mon. 
Not. R. Astr. Soc, 226 (1987) 739. 

3. D.W. Schwenke, J. Chem. Phys., 89 (1988) 
2076; 1990, J. Chem. Phys., 92, 7267 

4. D.R. Flower, Mon. Not. R. Astron. Soc, 297 
(1998) 334. 




500 1000 1500 2000 

T[K] 



500 1000 1500 2000 

T[K] 



Figure 4. Thermal rates calculated with the DJ 
PES for the 00 -> 22, 02 -> 22, 20 -> 00 and 
11 — > 13 transitions and other results |2I4| . 



5. D.R. Flower, E. Roueff, J. Phys. B: At. Mol. 
Opt. Phys., 31 (1998) 2935. 

6. S.K. Pogrebnya, D.C. Clary, Chem. Phys. 
Lett., 363 (2002) 523. 

7. S.Y. Lin, H. Guo, Chem. Phys. 289 (2003) 
191. 

8. M. Bartolomei, M.I. Hernandez, J. Campos- 
Martinez, J. Chem. Phys., 122 (2005) 064305. 

9. B. Mate, F. Thibault, G. Tejeda, J.M. Fer- 
nandez, S. Montero, J. Chem. Phys., 122 
(2005) 064313. 

10. F. Gatti, F. Otto, S. Sukiasyan, H.-D. Meyer, 
J. Chem. Phys. 123 (2005) 174311. 

11. R.A. Sultanov, D. Guster, Chemical Physics, 
in press 2006. 

12. P. Dicp, J.K. Johnson, J. Chem. Phys., 113 
(2000) 3480; ibid. 112 (2000) 4465. 

13. A.I. Boothroyd, P.G. Martin, W.J. Keogh, 
M.J. Peterson, J. Chem. Phys., 116 (2002) 
666. 

14. J. Le Bourlot, G. Pineau des Forets, D.R. 
Flower, Mon. Not. R. Astron. Soc, 305 (1999) 
802. 

15. J.M. Hutson, S. Green, MOLSCAT VER. 14 
(1994) (Distributed by Collabor. Comp. Proj. 
6, Daresbury Lab., UK, Eng. Phys. Sci. Res. 
Council, 1994) 



